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The surprising results by the BarBar collaboration on the ny transition form factor require new 
thoughts about the high-Q 2 dependence of the form factors with virtual photons. We make use 
of the anomaly sum rule [J. Horejsi and O. Teryaev, Z. Phys. C 65, 691 (1995).] which relates 
the hadron spectral density to the axial anomaly [S. Adler, Phys. Rev. 177, 2426 (1969); J. S. 
Bell and R. Jackiw, Nuovo Cimento A 60, 47 (1969) .]We study the quark-hadron duality relation 
for this sum rule and find out that the increase of the rescaled form factor Q 2 F- K1 (Q 2 ) ~ log(Q 2 ) 
suggested by the BaBar data requires the presence of a 1/ s-correction term in the relation between 
the one-loop spectral density and the hadron-continuum spectral density. 
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1. MOTIVATION AND RESULTS 

The processes 7*7 — > 7t(t7, 77') are of great interest for 
our understanding of QCD and the meson structure. In 
particular, the dependence on the space-like momentum 
Q of the incoming virtualphoton has been the subject 
of detailed experimental pHa] and extensive theoretical 
investigations (recent references see )■ QCD fac- 

torization predicts for large Q 2 the asymptotic behaviour 
17] Q 2 F^(Q 2 ) -+ V2fn, with / w = 0.131 GeV denoting 
the pion decay constant. Similar relations follow for 77 
and 7/ after taking particle mixing into account. Within 
errors, this saturation property is indeed observed for the 
1) and rj form factors. However, recent high-Q 2 data up 
to 35 GeV 2 indicate that Q 2 F V1 (Q 2 ) does not saturate at 
large Q 2 but increases further. FigO]compares these data 
with the theoretical formula obtained by Brodsky and 
Lepage [TtJ which interpolates the values of Ft T1 {Q 2 ) at 
Q 2 = given by the axial anomaly with the asymptotic 
form mentioned above. So far, no compelling theoretical 
explanation of the qualitatively different behaviour of the 
7T7 form factor compared to rjj and 7/7 form factors has 
been found; as concluded in 0, [HI, EH, 0] the behaviour 
of the 7T7 form factor is hard to explain in QCD. 

The exact anomaly sum rule [18| offers an interest- 
ing possibility [lpj to consider the transition form factor 
without directly referring to QCD factorization theorems. 
The aim of this letter is to address the Q 2 -behaviour of 
the P7 transition form factor from this perspective. 

Our starting point for a treatment of the three pro- 
cesses is the investigation of the amplitude 

(0|#7(«&)7*(«i)) = e% ia p(p\ qi ,q 2 ) e?4, 

p = qi + Q2- (l) 

Here E\ t 2 denote the photon polarization vectors. This 
amplitude is considered here for q\ = —Q 2 and q\ = 0. 
Its general decomposition contains four invariant form 
factors, but for our purpose only the single form factor 



with the Lorentz structure proportional to p^ is needed 

M 

Tp a p(p\qi,q2) = P^aHq l q 2 iF{p 2 ,Q 2 ) + ... (2) 

The invariant amplitude F(p 2 ,Q 2 ) may be written in 
terms of its spectral representations in p 2 at fixed Q 2 : 

F( P 2 ,Q 2 ) = - [ ^A( S ,Q 2 ). (3) 

7T J S — p Z 

4m 2 

In perturbation theory one obtains the spectral density 
as an expansion in powers of a s : 

A QC d(s,Q 2 ) = 

A^r>,Q 2 ) + ^A« D ( S ,Q 2 ) + O(o 2 ), (4) 

where the lowest order contribution Aqq D (s,Q 2 ) is ob- 
tained from the triangle diagram with the axial current 

Q 2 F ny {Q 2 ) [GeV] 




Fig. 1: Form factor F nj vs Q 2 : experimental data from [J, 0] 
(black dots) and [j] (red squares); dashed line represents the 
interpolation [13] which coincides with the local-duality model 
of pjl], the solid line represents our fit described below. 
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and two vector currents in the vertices 118-2 
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y/\ — 4m 2 / s. 



(5) 



Here m denotes the mass of the quark propagating in 
the loop. The integral of Aq^ d (s,(5 2 ) from s = Am 2 to 
infinity has the remarkable property to be independent 
of Q 2 and m [TJ[ and gives the axial anomaly (2l1 | 



dsA 



<°> ( s O 2 ) - J_ 

QCDl*' ^ ' _ 27r' 



(6) 
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The Adler-Bardeen theorem [22[ states that radiative cor- 
rections to the anomaly vanish, requiring that 



ds A 



(i) 

QCD 



(s,Q 2 



0, i > 1. 



(7) 



4m- 



This relation however does not require that all higher- 
order spectral densities Aq CD (s, Q 2 ) vanish; only the in- 
tegrals should be zero. The two-loop spectral density 
Aqq D (s,Q 2 ) was found to be identically zero [H, [24]]. 
Higher-order spectral densities have not been calculated. 
However, arguments have been given that all higher-order 
spectral densities may not vanish identically jl4j |. 

Non perturbative QCD interactions will strongly influ- 
ence F(p 2 , Q 2 ) producing a meson pole (and, correspond- 
ingly, a (near) delta function of its absorptive part) and 
a hadron continuum. Nevertheless, the integral over the 
entire absorptive part A(s,Q 2 ) remains unchanged. It 
still represents the anomaly: 



oo oo 

JdsA(s,Q 2 ) = J dsA QCD ( s ,Q 2 



1 

2tt' 



(8) 
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For the case of the isovector 



uu~dd 
V2 



axial current, the 

spectrum contains the 7r°-meson. Thus, the absorptive 
part of F(p 2 , Q 2 ) has the form 



A( S ,Q 2 ) 



(9) 



n6(s - ml) V2U F^{Q 2 ) + 6{s ~ s th ) A^(s, Q 2 ). 

Here A^ t (s,Q 2 ) denotes the hadron-continuum contri- 
bution in the isovector channel. Ft T1 {Q 2 ) then takes the 
form 



Fffj {Q ) 



2\/2 TT 2 /^ 



DC- 

l-tojds A^(s,Q 2 ) 



(10) 



For the 777 and 7/7 form factors, one needs to consider the 
isoscalar currents qq — (uu + dd)/y/2 and ss, separately: 



Fa q (Q 2 ) 



F S s{Q 2 ) 



1 



2\/2 7T 2 /„ 



2V2 7T 2 / S 
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1-2ttJ d s Al s oat (s,Q 2 ) 

Sth 



(11) 



For each channel, the relevant threshold Sth should be 
used. Taking 77 — 77' mixing (25l . l26j into account and 
using the corresponding quark charges one finds 



Frff {Q ) 
Fr)>-t(Q ) = 



-^=F qq (Q 2 ) cosc}>-^F Ss (Q 2 ) sin0, 

5 1 
-^-j=F qq {Q 2 ) sin^-f -F- SS {Q 2 ) cos</>. 



(12) 



The r\ — rj' mixing angle <f> is known to be ~ 39°; the 
decay constants are taken to be f q = 1.07 , f s = 1.36/^ 

M- 

According to (fT0|) and (fTT|) . the calculation of the P7 
form factors requires an Ansatz for the continuum spec- 
tral densities A cont (s,<3 2 ) for all three cases. 

At this place the duality concept may be used by re- 
placing the integrand from threshold to infinity by the in- 
tegrand obtained from perturbative QCD. For large val- 
ues of s, above the resonance region, the spectral density 
is anyhow expected to be very well represented by pertur- 
bative QCD. At lower s, however, the continuum spectral 
density is strongly distorted by strong interactions and 
may not be fully accounted for by the rapid onset of the 
perturbative expression. Moreover, the thresholds should 
be in general replaced by effective thresholds. This holds 
in particular for the qq and ss dispersion representations 
which determine the rj and rj form factors. Below the 
corresponding effective thresholds, the isoscalar interac- 
tion is very weak and is expected to be negligible. 

A simple Ansatz for A(s, Q 2 ) which can to some extent 
parametrize these effects is to write 



A cont (s, Q 2 ) = 9{s - s th )R{s)A ( V{s, Q 2 ), 



with R(s) 



1 



(13) 



Finite value of r, < r < s t h, smoothens the threshold 
behavior. Moreover, the r/s term strongly affects the 
high-Q 2 dependence of the form factor. We obtain this 
way a logarithmic increase with Q 2 of the form factor 
multiplied by Q 2 : 

Q 2 F(Q 2 ) ~ -&( Sth -r)+rlog ( 9l±I±) . (14) 



Sth 



Sth 



To derive (fT4"|) . we set m — (see the remark below). 
One may of course assume a more sophisticated behav- 
ior in the near-threshold region; in fact, the details of 
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Fig. 2: Form factors Fp 7 (P = r], rj') vs Q 2 : experimental data 
from P, 0| (black dots), [fj (red squares), and the data bor- 
rowed from the time- like region 0] (green triangle). Dashed 
line - results from local-duality model of [3, solid line - our 
fit. 



A con t(s, Q 2 ) at low s influence the form factor at small 
Q 2 . However, these details do not change the l/Q 2 - 
behaviour of Fp 7 (Q2) at large Q 2 . Our crucial observa- 
tion is that a wanted ~ log(Q 2 ) behavior of Q 2 Fp 1 (Q 2 ) 
requires the presence of the r/s term in R(s) (|13j) . 

Making use of the formula (TUT) , one finds values of the 
fit parameters Sth and r which lead to a good description 
of all existing experimental data, Figs. [T] and [5] These 
parameters are as follows: for the pion case we use the 
physical threshold value s t h = (Sm^) 2 = 0.165 GeV 2 and 
set also r = (3m„-) 2 ; then the continuum spectral density 
vanishes at sth- 

It is less obvious which parameter values for the r\ 
and rj' cases should be used because their structure is 
influenced by the strong gluon anomaly (26|. We set 
s t h = 0.56 GeV 2 for the qq channel and s t h = 0.76 GeV 2 
for the ss channel. These values are close to the {11} 
and the |22| elements of the rj — rj 1 mass matrix, respec- 
tively |26|. The fit for the r parameters yields r = 0.05 
GeV 2 for both channels. A good description of the data 
is achieved, see Fig. [2j Notice that this value of r is 
~ 1/3 of the corresponding value in the 1=1 channel. 
With regard to the quark mass values, one can safely set 
m = in (TT31 . at least for the it- meson case. For the 



Fq q and F Ss amplitudes needed for the rj and rf form fac- 
tors the dependence on possible values for m q and the 
strange mass value m s can easily be checked. As long 
as one sticks to current quark masses there is no visible 
difference within present accuracies from setting also for 
these cases m q = m s = 0. Therefore, the integrals in- 
volved above can be calculated analytically. The result 
has already been shown in (TH| . 

2. DISCUSSION 

We studied the 7*7 — > P transition form factors by 
means of the exact anomaly sum rule which relates the 
integral over the hadron spectrum to the axial anomaly. 
In order to isolate the pseudoscalar mesons from this sum 
rule, we make use of quark-hadron duality in the follow- 
ing way: we replace the integral over the continuum by 
the integral over the spectral density given by the di- 
agrams of perturbation theory. We include, however — 
essentially, by hand — an additional multiplication factor 
R{s) = (l — I). This factor goes to unity for high val- 
ues of s above the resonance region where perturbative 
QCD should fully describe the continuum. This multi- 
plication factor serves two purposes: (a) It simulates a 
better threshold behavior and (b) it modifies the high-Q 2 
dependence of the integral over the continuum in a way 
suggested by recent high-Q 2 measurements. Our main 
results are as follows: 

(i) The strong difference between the hadron spectral 
density and the perturbative QCD spectral density in 
the low-s region influences only the low-Q 2 structure, but 
not the large-Q 2 asymptotics of the form factor obtained 
from the anomaly sum rule. Using the parameter r dif- 
ferent from zero and positive, one easily obtains a good 
agreement for the low-Q 2 data, using threshold values 
well motivated by phenomenology. However, the increase 
observed for the 7r7 form factor by BaBar cannot be un- 
derstood by subtleties of the behaviour of the continuum 
spectral density near threshold. 

(ii) The deviation of the large-Q 2 behaviour of the form 
factor F 7Ty (Q 2 ) from the 1/Q 2 asymptotics predicted by 
QCD factorization requires a deviation of order 0(r/s) 
between the hadron-continuum spectral density and the 
lowest-order perturbative QCD spectral density in the 
region of larger s values. The multiplication factor R(s) 
in (fT3"|) of the general form 

R(s) = 1 - - + £i + . . . (15) 
s s z 

leads to the form factor which behaves at large Q 2 as 

Q 2 F P7 (Q 2 ) ~rlog(Q 2 / Sth ) (16) 
+ (sth - r - ri/ath) + O (log(Q 2 )/Q 2 ) . 

The logarithmic increase of Q 2 Fp 7 (Q 2 ) occurs propor- 
tional to r. Terms of higher order of l/s in R(s) do not 
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contribute to the leading log(Q 2 )/Q 2 -behaviour of the 
form factor. 

(iii) The BaBar data on Ft T1 (Q 2 ) are well described by 
using s t h = (3m„.) 2 , the true threshold for starting the 
continuum integral, and by setting r = (Sm^) 2 , a con- 
stant of the formal order of the light-quark mass. Since 
a good description of the data is obtained in this way, 
the theoretical treatment of the underlying (AW) cor- 
relator in the chiral limit appears to be insufficient. Fur- 
ther studies of the effects beyond the chiral limit may 
be needed. Moreover, one should better understand the 
possible origin of such 0(r/s) terms within QCD. So far 
it is a mere conjecture suggested by the data. 

(iv) Obviously, the proposed 0(r/s) term in R(s) which 
modifies the perturbative QCD spectral density leads to 



the violation of the QCD-factorization theorem. Only if 
the correction factor R(s) is replaced by unity at some 
high value of s, the form factors multiplied by Q 2 will 
saturate at Q 2 larger than this s-value. Since we can- 
not say where this will occur - if it occurs at all - the 
true asymptotic behavior for very large Q 2 cannot be 
predicted. 
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